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During the reporting period, we submitted 9 papers for publication and 
presented 1 review talk at a scientific meeting. I'he papers and review talk 
are listed below, and abstracts of the papers are attached. 

The papers are in the areas of plasma expansion phenomena, the polar 
wind, and auroral plasma processes. In one paper, we studied the development 
of electric fields in an expanding plasma (paper 4). These electric fields 
accelerate ions to suprathermal energies and are similar to those in double 
layers with a zero current. With regard to the polar wind, we showed that hot 
magnetospherlc electrons (polar rain, polar showers, and polar squall) have a 
pronounced effecf on the polar wind (paper 1). These hot electrons act to 
increase the outflow velocity and energy, with the energy gain being as 
large as 1 to 2 keV. We also showed that, contrary to previous pripdictions, 
there is no 0^ charge exchange barrier and substantial fluxes of O"*" ions can 
escape with the polar wind (paper 3). 

In the auroral plasma physics area, we studied the excitation of elec- 
trostatic waves by field-aligned auroral electron beams (paper 2). We also 
showed that the auroral field-aligned current density can be large enough to 
excite Buneman double layers (paper 5). For situations that lead to strong 
double layers, we showed that the temporal evolution of the potential profile 
is controlled by current fluctuations (paper 6). In addition, we conducted 
two-dimensional partlcle-ln-cell simulations and studied the high frequency 
wave turbulence excited by an auroral electron beam of finite width perpendic- 
ular to an ambient magnetic field (paper 7), and we studied the formation of 
V-shaped auroral potential structures (paper 8). Finally, we wrote a review 
article on numerical simulations of double layers and auroral electric fields 
(paper 9). 


Presentations 


R.W. Schunk, A Global Ionospheric Model, Invited Review , Presented at the 
"Optical Ground-Based Aeronoray Workshop II”, June 20-22, 1984; Ann 
Arbor, Michigan. 
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Abstract 

A seraikinetlc model was used to study the effect that hot electron popu- 
lations have on the polar wind. The model was used to describe the steady 

state plasma flow along diverging geomagnetic field lines in the collisionless 

•f* *4“ 

regime at high altitudes. The plasma contained 0 and U ions and both hot 
and cold electron populations. Several hot electron populations were consi- 
dered, including the polar rain, polar showers, and polar squall. Estimates 
of hot electron parameters based on characteristic energy and flux measure- 
ments Indicate that the hot/cold electron temperature ratio varies from 
10 - 10^ and that the percentage of hot electrons varies from 0.1% to 10% at 
4500 km. For ratios at the lower ends of these ranges, the polar wind solu- 
tions with hot electrons are similar to those obtained previously for super- 
sonic H outflow without hot electrons. For higher hot electron temperatures 
and a greater percentage of hot electrons, there is a discontinuity in the 
kinetic solution, which indicates the presence of a sharp transition. This 
traiVb'ition corresponds to a contact surface between the hot and cold elec- 
trons- Along this surface, a double layer potential barrier exists, which 
reflects the cold ionospheric electrons and prevents their escape. The pres- 

4 - 

ence of the hot electrons acts to increase the supersonic H outflow velocity 

+ + 

and H energy, but does not affect the already saturated 11 escape flux. The 

J- 

H energy gain may be as large as 1 to 2 keV . With regard to 0 , the hot 
electrons act to reduce the potential barrier, thereby allowing more 0 ions 
to escape. A significant enhancement in the escape flux can occur depend- 
ing on the hot electron density and temperature. 
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Abstract 

The excitation of waves at frequencies u < where Is the lower 
hybrid (AH) frequency, by a low energy electron beam Is examined under the 
electrostatic approximation; the frequency range w < includes the 
ion-cyclotron frequencies (n^^) and their, harmonics in a multi-ion plasma. 
Since the beam velocity » v^^, the electron thermal velocity, the excited 
waves have phase velocities Thus, the background electrons 
can be treated as cold and the electron density fluctuations are small. In 
the limit of cold ' ion approximation, the excited waves are either the 
ion-cyclotron resonance cone (ICRC) waves or the beam-plasma (b-P) mode. Both 
of these waves show stop bands near the gyrofrequencies of the constituent 
ions. The ICRC waves are limited to the highest gyrofrequency while the B-P 
mode propagates even beyond it. A comparison of the characteristics of the 
electrostatic ICRC waves with those determined from a full electromagnetic 
(FEM) treatment at frequencies u < dhows that the electrostatic approxima- 

-ii-i. '.Mlid 'when beam energy is low (-lOOeV) and for relatively more ener- 
getic beams the validity holds good only when the perpendicular wave number 
kj^ is sufficiently large (k_L>10 n(H )/C, n being the H gyrofrequency, and C 
the velocity of light in vacuum). Warm ion effects on the waves are examined; 
when ^ ^ 0.3, where r^(H'^) is the h"^ Larmor radius in a He^-H^ 
plasma, the excited waves are the EIG waves. Because of their phase velocity 

V » V , these EIC waves are pure ion-Bernstein (PIB) waves. The charac- 
p te 

teristics of the electron-beam excited PIB waves are compared with their coun- 
terpart, the neutralized ion-Bernstein (NIB) waves, which are EIC waves excit- 
ed by drifting bulk electrons. The consideration of electrostatic waves in 
the frequency band 0 (h"^) < u < u shows that the B-P waves with 



« 1 have considerably larger temporal growth rates than the EIC 
waves. The interconnections between LH and B-P waves are found to be similar 
to those between EIC and B-P waves; B-P waves have long perpendicular wavel- 
engths (Xj_ ) while EIC nnd LH waves have comparatively much shorter ones. The 
relevance of ICRC and PIB waves to satellite observations of waves in the ELF 
band is discussed. The B-P waves appear to be similar to the electrostatic 
broadband noise (BEN) emissions in the frequency range from below the 
ion-cyclotron frequency to the lower hybrid frequency observed from various 

' * V 

satellites in the auroral and cusp regions. 
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Abstract 

The continuity and momentum equations for ionospheric 0 were solved for 

typical polar cap condition's. Density profiles were, obtained both with and 

4 * 

without allowance for accldentally-resonant charge exchange (ARCS) between 0 

4 * 

and H. It was found that ARCE acts to reduce the maximum (or limiting) 0 
escape flux by less than 30% for typical atmospheric conditions (T > 1000 

03 M 

K) and by only a factor of about 3.5 for atmospheres with low exospheric tem- 
peratures (T^ 750 K) . With allowance for ARCE the limiting 0*^ escape flux 

is of the order of 5 x 10 to 10 cm s , depending on the atmospheric con- 
ditions, which easily accounts for the existing measurements. This indicates 
that, contrary to previous predictions, there is no 0^ charge exchange bar- 
rier, and consequently, it is not necessary to have an acceleration mechanism 

•f* 

at low altitudes in order to have an appreciable 0 escape flux. Instead, an 
4 * 

0 depletion at high altitudes due to some acceleration mechanism is suffi- 

4* 

dent to cause large 0 fluxes to escape from the ionosphere to the magnetos- 
phere . 
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The phenomenon of plasma expansion occurs in a variety of situations 
ranging from laser fusion t6 space physics. It is now believed that plasma 
expansion plays an important role in the electrodynamics of the wake of bodies 
in plasmas. In space plasmas, these bodies can be either man-made or natural 
objects, such as the terrestrial moon and planets, etc. The phenomenon of 
plasma expansion has been fouiid to be relevant to the initial expansion of the 
polar wind and to Interhemispherlc plasma flows. An important characteristic 
of a plasma expansion is the development of polarization electric fields, 
which can be localized, 'iiivise electric fields accelerate ions to suprathermal 
energies. Since the electrons and ions in the expansion drift with equal vel- 
ocities, there is no net current. Thus, the localized electric fields asso- 
ciated with expanding plasmas are like double layers with a zero current. 
These double layers move with the expanding plasma front- The characteristics 
of the localized electric fields in both single and multi-ion plasma expan- 
sions, as seen from numerical simulations, will be presented. These charac- 
teristics will be compared with those of classical double layers with cur- 
rents. 
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Results from both one'-dlmcnslonal numerical simulations and space mcas-* 
urements on the magnitude of the field-aligned currents will be presented to 
show that Buneman double layers can be driven on atroral field lines. The 
temporal evolution of the current density and electric field will bo used to 
show the various stages leading to the formation of double layers. The cur- 
rent density in the plasma displays .a cyclic interruption and recovery. The 
Buneman instability is triggered during current recovery. The current interr- 
uption occurs along with the double layer formation, which is also cyclic. 
Some dynamical aspects of Buneman double layers i^ll be discussed. It will be 
shown that when the current through the plasma is relatively large, a common 
feature of the turbulence on the high potential side of the double layers is 
the formation of electron holes. On the other hand, when the current drops 
below the thermal current, a localized potential dip foians on the low poten- 
tial end of the Buneman double layers. 


CURRENT FLUCTUATIONS AND DYNAMICAL FEATURES OF 
DOUBLE UYER POTENTIAL PROFILES 

Nt Singh, H. Thlcnnnn and Ri W. Schunk 

Center for Atmospheric and Space Sciences . 

Utah State University, Logan, Utah 84322 USA 

I 

The low frequency oscillations caused by potential relaxation instabili- 
ties (PRI) occurring in one-dimensional simulations of strong double layers 
and the accompanying formation of potential dips on the low potential end of 
the double layers is reviewed. T\JO-dlmensional (2-d) numerical simulations of 
double layers in a magnetized plasma show oscillations in the current caused 
by both PRI and electrostatic ion cyclotron (EIC) waves. During current 
troughs, the double layers tend to develop localized dips at the low potential 
end of the parallel profiles, and their strength tends to increase. The 
occurrence of potential dips in association with current Interruptions in very 
diverse situations shows that such a feature is universal. 
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High frequency (u)~ti>pe< Gg ) wave turbulence excited by an 
electron beam of finite width perpendicular to an ambient 
magnetic field in a background plasma is described. The 
results discussed here are from a 2 1/2-dimensional numeri- 
cal simulation with a plasma bounded in the direction of the 
current flow and unbounded with the periodic boundary condi- . 
tion in a direction transverse to it. During the early stages 
of the beam-plasma interaction \-a very intense turbulence occurs, 
which eventually disrupts and at a late stage a quasi-steady 
state with an appreciable level of turbulence follows. The high 
frequency turbulences undergo considerable modulations by va- 
ridus low frequency phenomena. Although the . electron 
current £s limited to a narrow channel’, . the v/avea excited 
by the current sprea.ds efficiently outside the channel. 
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SOKE FEATURES OF AURORAL ELECTRIC FIELDS AS SEEN 
IN 2D numerical SIMULATIONS 
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ABSTRACT 

Reoults of 2D plaeno elmulatlona, which ahow the formation of V-ahaped double laycra with 
larger perpendicular than parallel electric field amplltudce, the excitation of clectroatotlc 
Ion-cyclotron wavea and the occurrence of return currento, arc preoented and related to auro- 
ral obaervatlona. 
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ABSTRACT . 

Recent niaierieal alBulatlona on double lapera XOU) are reviewed and dlecuaaed In the light of 
relevant apace obaervatlona In the auroral plaana< Tvo-dlnenilonal (2-d) DLe driven bjr cur- 
rent aheeti of a finite thlckncaa L ahov different charaeterlatlca depending on whether l' < 
P( or t » pj, p, being the ion gyroradlua. In the fonaer eaae, V-ahaped DLa fora with near^ 
if e<}ual parallel and perpendicular potential dropa { , which are dctemlncd by the teapera- 
turea of the varloua conatlcuent Xonai In the lacear'’caae, Che aajor parallel potential drop 
occura outalde the current aheec and the pctpendfpnlar elceCrlc flelda arc locallted at the 
edgea of Che current aheeCi Soae laporCant fcaturea, auch aa the aaplttudea and acale 
lengtha of Che parallel and perpendicular electric flelda, aeen In cheat alaulaClona are 
qualitatively alnllar to thoae obaerved In apace. The review of one-dlnenalonal aloulaclona 
covtra the followlngt (1) Ttie dynanlcal fcaturea of DLa and their potalble appllcaclona to 
fllckcrlog aurorae and low frequency fluctuaclona of Inverced-V eventa; (11) lon-acouaclc 
double layera with auatalned currenCa and their relatlonahlp to weak DLa obaerved Ij auroral 
plaanaa; (111) Recent DE 'obaervatlona of large field-aligned currenta (> 100 pA/* ) and the 
excitation of Buncaan double layera. 
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